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A time periodic driving on a topologically trivial system induces edge modes and topological 
properties. In this work we consider triplet and singlet superconductors subject to periodic variations 
of the chemical potential, spin-orbit coupling and magnetization, in both topologically trivial and 
nontrivial phases, and study their influence on the charge and spin currents that propagate along the 
edges of the two-dimensional system, for moderate to large driving frequencies. Currents associated 
with the edge modes are induced in the trivial phases and enhanced in the topological phases. In 
some cases there is a sign reversal of the currents as a consequence of the periodic driving. The edge 
states associated with the finite quasi-energy states at the edge of the Floquet zone are in general 
robust, while the stability of the zero quasi-energy states depends on the parameters. Also, the spin 
polarization of the Floquet spectrum quasi-energies is strong as for the unperturbed topological 
phases. It is found that in some cases the unperturbed edge states are immersed in a continuum of 
states due to the perturbation, particularly if the driving frequency is not large enough. However, 
their contribution to the edge currents and spin polarization is still significant. 

PACS numbers: 74.40.Gh, 72.25.-b, 74.40.Kb 


I. INTRODUCTION 

Topological systems have attracted great 
intereslP^l and, specifically in the case of topolog¬ 
ical superconductor^, a great effort has been put 
towards prediction and detection of Majorana fermion^. 
The robustness of topological systems to perturbations 
is a key feature with wide potential applications, 
such as the edge state Majorana fermions as building 
blocks for topological quantum computation, due to 
their non-abelian statistics. It has been shown that 
topological systems are quite robust to a quantum 
quench, as exemplified by the toric code model^. It 
has also been shown recently that, in an infinite system, 
the Chern number can not be changed by a unitary 
evolution, even though it is possible to change the Bott 
index and topology if the system has a finite siz^, as 
obtained earlier in the case of quenches in topological 
superconductor^ Therefore, the response to a time 
dependent perturbation of the edge states in finite sys¬ 
tems may not be protected by topology. Furthermore, 
quenches in superconducting systems with topological 
properties, performed self-consistentlji^, showed the 
importance of the topological properties in the evolution 
of the systerrPS! amj raised questions reg arding the 
survival of the topological order to the quenchPEi®^. So 
the issue is not resolved and is attracting considerable 
attention. 

Several examples that host Majorana fermions as edge 
states are p rovided by several superconducting p-wave 
systemsHUSl and various other systems that mix super¬ 
conducting order (eventually by proximity effects) with 
Zeeman fields and/or spin-orbit couplingISI or tunneling 
of Cooper pair^I^. A recent proposal for the existence of 
localized Majorana edge states at the ends of a magnetic 
chain in contact with a conventional superconductoili^ 


has been confirmed experimentalljff^ (there is however 
some debat e because of their unexpected very localized 
natur^2^121J) and it has been shown, combining topolog¬ 
ical with non-topological regions in wires, that it also 
provides a good basis of states to implement the braiding 
properties of the Majorana fermion^^, with potential in¬ 
terest in quantum computation. The fermion exchange in 
a unidimensional system recently proposecP^ is also pos¬ 
sible to implement with a magnetic chain. Various other 
proposals have been presented for the existence of Ma¬ 
jorana zero energy states including recently multiband 
system^^miSI^ 

where the non-trivial interband coupling 
leads to interesting topological properties, as shown also 
in the context of topological Kondo insulator^^ZHH 

In higher dimensions the edge modes are in general 
propagating in the directionss parallel to the edge and 
localized in the perpendicular direction. Their nature 
as protected states suggests the existence of spontaneous 
currents, either charge or spin currents, that have aroused 
considerable interest in several contexts such as in the 
context of spintronics, as a means to generate and manip¬ 
ulate spin currents. The anomalous Hall effect due to the 
presence of spin-orbit coupling has been considerecP^**^ 
as well a s the appearance of currents due to magnetic 
impuritie^SSEI] study of the currents generated at 

interfaces between topologica l insu lators and a super- 

condutor has been considerecP^Ml ^jjg charge 

and spin conductances. A manipulation using electrical 
means has been proposecPSl as well as their spin correla¬ 
tions using a gate potentiaP^. Chiral p—wave supercon¬ 
ductors are expected to show spontaneous edge currents. 
Some other chiral systems were, however, shown to not 
have spontaneous current^^^^^ but, with some admix¬ 
tures of different pairings, these arise together with some 
non-trivial spin polarizatioiP^. Even though the topic 
has been studied for some time, the experimental verifi- 
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cation of the spontaneous generation of edge currents is 
still a topic of considerable attention. It is expected that 
a strong candidate in the class of chiral p-wave supercon¬ 
ductors Sr2RuOi will exhibit spontaneous currentj-' ^^ l ^^ l 
Even though the material Sr2RuOi has been shown to 
break time reversal symmetry due to the existen ce of 
the Kerr effect!^, no edge surrents have been founcP^EH 
even though evidence for the edge states has been seen 
by in-plane tunneling spectroscopj®!. Various propos¬ 
als have been presented to explain the absence of the 
edge c urrent s^ using the result of an effective high Chern 
number®lSlor effects of disorder!^. Also, it has been pro¬ 
posed that the absence of the edge currents may be due 
to a spin singlet pairiM instead of the generally assumed 
p-wave triplet pairing^. 

As mentioned above, it is interesting to study the ro¬ 
bustness of the edge states to time dependent pertur¬ 
bations. In the context of the Creutz ladder, it was 
shown that the presence of edge states modifies the pro¬ 
cess of defect production expected from the Kibble-Zurek 
mechanism, leading in this problem to a scaling with the 
change rate with a non-universal critical exponent!^. A 
similar result was obtained for the one-dimensional su¬ 
perconducting Kitaev model, where it was shown that, 
although bulk states follow the Kibble-Zurek scaling, 
the produced defects for an edge state quench are quite 
anomalous and independent of the quench rat c i ^ l ^^l The 
behavior of edge states under an abrupt quantum quench 
has also been considered very recently in the context of 
a two-dimensional topological insulatoil^, where it was 
found that, in the sudden transition from the topological 
insulator to the trivial insulator phase, there is a collapse 
and revival of the edge states^. Similar results were 
obtained for the one-dimensional Kitaev modeP^, and 
also studying the signature of the Majoranas in the en¬ 
tanglement spectrunP^. Their dynamical formation and 
manipulation has been considered irP^ ancP^l. The ef¬ 
fect of a sudden quench of the parameters of the Hamil¬ 
tonian of a two-dimensional triplet superconductor has 
also been studied, and the robustness of the edge states 
was considerecP. In general it was found that the edge 
states decay due to the quench, even though in some 
cases they are quite robust, such as in the case of weak 
spin-orbit coupling, or when there are matching momen¬ 
tum states in the initial and final states. The effect of 
parity blocking on the dynamics of the edge modes has 
been considered recently in which case the dynamics is 
restricted if there is a change in fermion parity accross 
the quenctP^. On the other hand, the Majorana zero 
modes lead to some universal non-equilibrium signature 
in the Loschmidt echo with an universal exponent as¬ 
sociated with the algebraic decaj^^. Also, the dynam¬ 
ics of the tunneling into non-equilibrium edge states has 
been proposed as a possible signature of the existence of 
these state^^. Non-equilibrium situations also may allow 
the transport of Majorana edges states using extended 
gapless regions with a small but finite overlap with the 
MajoranasSl. Their effect has also been considered irP^. 


While quenches, either abrupt or slow, in general 
destabilize the edge states, topological phases can be in¬ 
duced by periodically driving the Hamiltonian of a non- 
topological s^em, such as shown before in topologi¬ 
cal insulatord^SHSSl jjj topological superconductors, 
with the appearance of Majorana fermion^SSHsl Their 
appearance in a one-dimensional p-wave superconductor 
was studied in Ref.^^ and in Ref.l^ introducing exter¬ 
nal periodic perturbations; the case of intrinsic periodic 
modulation was also considerecP^. The periodic driv¬ 
ing leads to new topological stated, and to a general¬ 
ization of the bulk-e dge correspondence, that reveals a 
richer struci^r^^UZil as compared with the equilibrium 
situatioiP^^^. Similarly, in topological superconductors 
new phases may be induced an d manip ulated due to the 
presence of the periodic drivin^^SEUZH 

In general the problem is complex due to problem of 
dissipation through coupling to baths. Their effect have 
recently been considered and their detection has been 
proposed using transport propertiePsHUl^ as well as mag¬ 
netization signatures using the magnetic fields the cur¬ 
rents should produc^^. 

In this work we compare the charge and spin edge cur¬ 
rents of the unperturbed Hamiltonian with those of the 
perturbed Hamiltonian for triplet superconductors. We 
consider a system with no coupling to its environment 
except for the coupling to the periodic perturbation. To 
simplify, we consider perturbations where only one of the 
parameters of the Hamiltonian is periodicaly changed, 
such as the hopping, spin-orbit coupling, chemical po¬ 
tential or magnetization. We mainly focus on the last 
two parameter changes since are easier to implement ex¬ 
perimentally. In section H we review the Hamiltonian 
considered and in section HI we review the theory of the 
Floquet states that result from a periodic perturbation. 
In section IV we illustrate the method by presenting so¬ 
lutions of the quasi-energies for different cases, as a func¬ 
tion of the perturbation amplitude and frequency, and 
verify the convergence of the truncation procedure of the 
Hamiltonian matrix of the Floquet problem. In section 
V we present results for the edge currents and the spin 
polarization of the states of system. We conclude with 
section VI. 


II. UNPERTURBED TOPOLOGICAL 
SUPERCONDUCTOR 

We consider two-dimensional superconductors that 
display topological phases. Standard examples are triplet 
superconductors with p-wave symmetry. We also con¬ 
sider the presence of spin-orbit coupling and a Zeeman 
term (magnetization that may be due to proximity ef¬ 
fect or intrinsic). In the case of a non-centrosymmetric 
system, since parity is no longer conserved, a spin sin¬ 
glet pairing is also possible. Therefore we also consider 
an admixture of a conventional s-wave component. Even 
without a triplet pairing component a conventional su- 
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perconductor may also display topological properties in 
the presence of magnetization. 

The model considered here was s tudied in various ref- 
serences before such as in Refs.li^*^. We write the Hamil¬ 
tonian for the bulk system as 


k 


ffo(fe) A(fe) 
At(fe) 



where V'l;, V'-fet-V’- h) and 


Hq = ekO-Q - + Hr . (2) 


Here, = —2t(coskx + cosky) — ^ is the kinetic part, 
t denotes the hopping parameter set in the following as 
the energy scale (< = 1), /i is the chemical potential, 
fc is a wave vector in the xy plane, and we have taken 
the lattice constant to be unity. Furthermore, is the 
Zeeman splitting term responsible for the magnetization, 
in t units. The Rashba spin-orbit term is written as 


Hr = s ■ a = a (sin kya^ — sin k^cfy) , (3) 


where a is measured in the same units. The matrices 
(Jx^Uy, Uz are the Pauli matrices acting on the spin sector, 
and (Jo is the 2x2 identity. The pairing matrix reads^ 


K = i(d ■ cr) Uy = 


da; idy dz -j- Ag 
dz Ac dx. id^i 


( 4 ) 


where A^ is the symmetric part and the vector d 
parametrizes the anti-symmetric part. 

The energy eigenvalues and eigenfunction may be ob¬ 
tained solving the Bogoliubov-de Gennes equations 


f Ho{k) A(fe) 

{A^ik) -HU-k) 


— ^k,n 


( Un{kA) 
'i') 

'^^n( T) 

V Vn{-k,l) 


The 4-component spinor can be written as 


( 5 ) 


( 6 ) 


We consider either a superconductor with d = 

d{siiiky, — sinkx,0) or d — dz{0,0,sinkx — isinky) with 

or without a contribution from the local s-wave pairing 
Ag. The first case applies if the spin-orbit is strong. 
In this case the pairing is alignecP^ along the spin-orbit 
vector s. This case is denoted by strong coupling case. 
Relaxing this restriction allows that the two vectors are 
not aligned. This case is denoted by weak spin-orbit 
coupling and has been considered before in the con¬ 
text of the anomalou s Hal l effect and the calculation of 
the Hall conductanc^^nEU strong-coupling case 

d = {dx,dy,dz) = {d/a)s. The second case, usually 
called p + ip superconductor, has been the focus of great 
attention since it has spontaneous chiral states due to its 



FIG. 1: (Color online) Phase diagram for the spin triplet 
superconductor with d = {dx,dy,0). The phase diagramfor a 
s-wave superconductor is very similar except in the vicinity 
of zero magnetization (limited approximately by the dotted 
lines) where the system is topologically trivial. In the case 
of a triplet superconductor with d = (OjO,^^) with a com¬ 
plex ^-component, the phase diagram is also similar with the 
exception of a region around zero magnetic field (limited by 
dashed lines) that is topological, with Chern number 2. 


time-reversal breaking structure. In contrast, in the ab¬ 
sence of the Zeeman term the first case does not break 
time-reversal symmetry (TRS) and in the topological 
phases the edge states have a helicoidal structure. The 
system then belongs to the symmetry class Dili where 
the topological invariant is a Z 2 indejP^. If the Zeeman 
term is finite or in the case of the p + ip pairing, TRS is 
broken and the system belongs to the symmetry class D. 
The topological invariant that characterizes this phase is 
the first Chern number C, and the system is said to be a 
Z topological superconductor. 

The Chern number is obtainable as a n int egral of the 
Berry curvature over the Brillouin zon d^^ l ^^ . Summing 
over the occup ied bands the Chern number has been 
calculatecpSEl] ji; jg convenient to calculate the Chern 
number by computing the flux of the Berry curvature 
over plaquetes in the Brillouin zon^^. 

The results in the parameter space for the various pair¬ 
ing symmetries are shown in Fig. [fusing the parameters 
a = 0.6, d = 0.6 and in the case of strong spin-orbit cou¬ 
pling. In this case, changing these parameters only leads 
to a quantitative change of the shape of the energy bands 
and no qualitative changes are observed in the topolog¬ 
ical properties. As shown befor^i^, the spin-orbit cou¬ 
pling does not change the topology. The phase diagram 
for the case of a s-wave superconductor is very similar 
except for values of the magnetization smaller than the 
amplitude of Ag. In these regimes the Chern number 
vanishes and the system is topologically trivial for all 
values of the chemical potential. Something similar hap¬ 
pens in the case of the p + ip superconductor in the sense 
that some regimes at values of the magnetization smaller 
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than the pairing amplitude, the Chern number is dif¬ 
ferent from the case of the dx , dy pairing. In these regimes 
it is finite, (7 = 2, due to the broken TRS, even in the 
absence of a finite magnetization. However, for some val¬ 
ues of the chemical potential, a finite magnetization will 
change the Chern number to zero (even though there are 
unprotected edge states). 

Due to the bulk-edge correspondence if the system is 
placed in a strip geometry and the system is in a topo¬ 
logically non-trivial phase, there are robust edge states, 
in a number of pairs given by the Chern number, if time 
reversal symmetry is broken. There are also counter- 
propagating edge states in the Z 2 phases even though 
the Chern number vanishes, as in the spin Hall effect. 
In these phases time reversal symmetry is preserved and 
the Kramers pairs of edge states give opposite contribu¬ 
tions to the Chern number. Interestingly, turning on the 
magnetization (Zeeman field) time reversal symmetry is 
broken and the edge states are no longer topologically 
protected. However, it was found that, even in regimes 
where (7 = 0, there are edge states, reminiscent of the 
edge states of the Z 2 phases, as long as the chemical 
potential satisfies |yLt| < At. The regimes where the mag¬ 
netization is zero and \fj.\ > At are topologically trivial 
and in a strip geometry have no edge states. 

III. PERIODIC DRIVING 


Here Hd(k) is of the form of the unperturbed Hamilto¬ 
nian but with only one non-vanishing term. Looking for 
a solution of the type 

V'(fc,t) =e-*"('=)‘$(/c,t) (9) 

and using that d>(fc, t) = $(fc, t+T), where T is the period 
(w = ‘2,'kIT) and /(wt) = /(w(t -I- T)), one gets that 

{k, t) - <^(k, t) = e{k)<P{k, t) (10) 

The time-independent quasi-energies e{k) are the eigen¬ 
values of the operator H{k,t) — and the function 
$(A:,t) the eigenfunction. Since this function is periodic, 
we can expand it as 

$(fc,t)=^<^„(fc)e*™‘"‘ (11) 

ra 

Inserting this expansion in equation we obtain the eigen- 
system 

'^Hmm'ik)4>m'{k) = e{k)4>m{k) (12) 


The time evolution of a state under the influence of a 
time dependent Hamiltonian is given by 

t) = H{k, t)V'(fc, t) (7) 

where k is the momentum, t the time and h = 1. We can 
decompose the Hamiltonian in two terms a time indepen¬ 
dent one, H{k), and an extra term due to the external 
time-dependent perturbation, that we want to take as 
periodic with a given frequency, w, 

H{k,t) = H{k) + f{ujt)Hd(k) (8) 


The Hamiltonian matrix is given by 


Hmm'{k) = dmm'rnw + — dtc 


*H{k, t)e 


(13) 


Choosing a perturbation of the type /(wt) = cos(wt) 
the second term of the Hamiltonian matrix reduces to 
1/2 {Sjri' + l,m A~ — 


The time evolution of the state is then obtained solv¬ 
ing for the quasi-energies, e(k), and the functions 4>Tn{k) 
diagonalizing the infinite matrix 


/ 


V 


{m — 2)uj + H(k) ^Hd{k) 0 

^Hdik) {m-l)uj + H{k) \Hd{k) 

0 lHd{k) muj + H{k) 

0 0 ^^Hd{k) 

0 0 0 


0 

0 

\Hdik) 

(to -I- 1 ) 0 ; -I- H{k) 
\Hd{k) 


(to -I- 


0 

0 

0 

\Hd{k) 

2)uj + H{k) 


\ 


(14) 


/ 


The matrix can be reduced if the frequency is high 
enough and only a few values of to are needed. The 
photon-dressing of the band structures due to the mixing 


of the bands is the important effect we consideil^. We 
will neglect here any photon emission/absorption pro¬ 
cesses that affect the occupation numbers of the elec- 
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irons. Equivalently one may consider that the photons 
are off-resonance. The effects of creation/absorption of 
photons if the driving frequency is not large enough were 
discussed befor^^. The matrix Hd{k) will be chosen 
as the unperturbed Hamiltonian, where only one of the 
parameters, hopping, chemical potential, spin-orbit cou¬ 
pling or magnetization will be considered to vary with 
time. The first three parameters preserve time reversal 
symmetry while the magnetization naturally breaks time 
reversal symmetry if the unperturbed Hamiltonian is in a 
regime with vanishing magnetization. Emphasis will be 
placed on the effects of varying the chemical potential or 
the magnetization which are easilly tuned externally. In 
this last case it has been determined befor^^ that even 
though the low energy states have a very low energy, they 
may not be strictly Majorana fermions since the eigen¬ 
values of the Floquet operator (time evolution operator 
over one time period) are not strictly ± 1 . 

Due to the periodicity of the eigenfunctions, -|- 

T) = $(fc,t), the action of the evolution operator, U{t), 
on a state over a period, T, leads to the same state minus 
a phase 

mT)) =lt{T)\m) = (15) 

Therefore, the quasi-energies are defined minus a shift of 
a multiple oi w = 2TTjT, and we can restrict the quasi¬ 
energies to the first Floquet zone, defined by the interval 
—wj2 < e < w/2. States with quasi-energies e = w/2 
and e = —w /2 are therefore equivalent and there is a 
reflection of any bands as one exits the Floquet zone 
from above (or below) and as one enters from below (or 
above). Considering the particle-hole symmetry of a su¬ 
perconductor, 7 _g = 7 ! and the equivalence between the 
energies e = —w/2, w/2 one expects a new type of Ma¬ 
jorana mode in addition to any zero energy states, the 
usual Majorana modes. 

For the large enough driving frequencies considered 
here, the frequency description considered is particularly 
convenient. In the regime of small driving frequencies a 
time description is more convenient. 

IV. ENERGY BANDS AND QUASI-ENERGIES 

The energy spectra of the unperturbed superconduc¬ 
tors both in the trivial and in the topological phases have 
been determined before. The solutions for the quasi¬ 
energies of the perturbed Hamiltonians lead to bands 
that have a similar structure. Here we focus on the en¬ 
ergy and quasi-energy spectra in a ribbon geometry to 
give evidence to the states along the edges of the system. 

At large frequencies, w > At, the size of the trun¬ 
cated matrix is relatively small and the quasi-energies 
and physical properties (calculated over the first Floquet 
zone) converge fast for small values of m. Considering 
m = 0 one reproduces the Hamiltonian of the unper¬ 
turbed superconductor. The first approximation for the 
driven system is obtained considering to = 1 , 0, —1, then 


TO = 2,1,0,-1,-2 and so on. One may therefore use 
a short notation for the number of terms considered in 
the diagonalization of the Hamiltonian matrix by using 
I = 0,1,2,3, •••. The unperturbed case is denoted by 
I = 0 and the perturbed cases by I = 1,2, • • • (consider¬ 
ing that we are using 21 + 1 states). If the frequency w 
is small, one needs to consider large values of I and the 
problem of finding the edge states in a ribbon geometry 
quickly becomes heavy computationally. Increasing the 
value of the frequency it is easy to find that it is enough 
to consider I = 2, since taking I = 3 leads to very similar 
results, with a good accuracy. 

To illustrate the structure of the quasi-energies we con¬ 
sider first the case of / = 1 and a dx,dy superconductor 
in a topologically trivial phase. The results are shown 
in Figs. and where we consider periodic drivings 
in the hopping and in the magnetization for moderate 
couplings of id = Mzd = 2. We consider frequencies 
w = 4,6,8,12. As the frequency increases edge states 
appear at the border of the Floquet zone around w/2 
(and —w/2). As the frequency increases further there is 
a continuum of states that fills the gap around zero en¬ 
ergy. On the other hand, at higher energies close to the 
Floquet zone edge, as the frequency increases the edge 
states appear on top of the continuum and prevail as the 
gap around w/2 opens. Therefore, as a consequence of 
the periodic driving edge states appear due to the Flo¬ 
quet spectrum in a system that is trivial if unperturbed. 
The edge states around w/2 visible for the higher fre¬ 
quencies are robust and converge as I increases. The two 
type of perturbations couple to states differently partic¬ 
ularly in the spin structure. Also, the hopping preserves 
time reversal symmetry while the magnetization breaks 
time reversal symmetry. However, their effects at this or¬ 
der are very similar with respect to the edge states. The 
main difference is the response of the bulk states to the 
two types of perturbations, particularly seen at moderate 
frequencies where in the case of the magnetization pertur¬ 
bation the continuum appears at the energies where the 
edge states appear, even though at different momenta. 

In Fig. we show the convergence of the truncation 
procedure of the Hamiltonian matrix for a high value of 
the frequency w = 6. We show results for I = 1, 2, 3 and 
see that the results for I = 2 and I = 3 are very similar. 
For a smaller frequency value, w = 4, the convergence 
is also fast. Moreover, for the set of parameters consid¬ 
ered, with a relatively high coupling of M^d = 4 we see 
that, in addition to a fast convergence and effectiveness 
of the truncation procedure, besides the edge Majorana 
states close to the Floquet boundary robust Majoranas 
also appear at zero energy. For smaller frequencies and 
weak couplings the convergence is in general slower. As 
a consequence, and to reduce the computational effort, 
we consider a frequency w = 6 in most results and limit 
the results to I = 2. 

It is also interesting to determine the effect of the peri¬ 
odic driving when the system is originally in a nontrivial 
topological phase. In Fig. [^we show the effect of driving 
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FIG. 2: Floquet spectra for Z = 1 for a dx,dy triplet superconductor in a topologically trivial phase with a = 0, Mz = 0, /i = 
— 5, As = 0, dz = Q,d — 0.6 where the hopping is changed with time with different frequencies w = 4,6,8,12. The periodic 
driving is tdcoswt with td = 2. Note that the spectrum is symmetric around zero quasi-energy. 



FIG. 3: Floquet spectra for Z = 1 for a dx,dy triplet superconductor in a topologically trivial phase with a = 0, Mz = 0, /i = 
— 5, As = Ojdz = 0, d = 0.6 where the magnetization is changed with time with different frequencies w = 4,6,8,12. The 
periodic driving is Mzd cos wt with Mzd = 2. 
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FIG. 4: Evolution of Floquet spectra for various values of I, I — 0,1, 2, 3 for a. dx,dy triplet superconductor in a topologically 
trivial phase a = 0.1, = 0, /r = —5, As = 0.1, = 0,d = 0.6 where the magnetization is changed with time with frequency 

w = 6. The periodic driving is M^d cos wt with M^d = 4. The case Z = 0 is the unperturbed Hamiltonian. The results show 
the stability of two sets of Majorana modes both at zero energy and at the limit of the Floquet zone. 



FIG. 5: Evolution of Floquet spectra for I — 1 for a d^,dy triplet superconductor in a topologically non-trivial phase 
a = 0,Mz = 2,/i = —5, As = 0,dz — 0,d = 0.6 where the hopping is changed with time with frequencies w = 6,12. The 
periodic driving is td cosuit with td ~ 0.5 (top row) and id = 2 (bottom row). 




























































































the hopping with different frequencies w = 6, 12. The re¬ 
sults are shown for Z = 1. We find that the continuum fills 
completely the gap, as illustrated in the case of rc = 6. 
Also, the same happens at the edge of the Floquet zone, 
particularly if the coupling is small = 0.5. Increas¬ 
ing the coupling to id = 2 an edge state emerges at the 
boundary of the Floquet zone. Also both for moderate 
{w = 6) and high frequencies {w = 12) edge states ap¬ 
pear near zero energy. However, the influence of the bulk 
states is considerable and at moderate frequency w = 6 
they are superimposed on the edge states. 

The topology of the system in the presence of the pe¬ 
riodic driving may also be considered in terms of the 
system in a torus. As mentioned above, the classification 
is richer and the bulk-edge correspondance is more sub¬ 
tle. It has been proposed^ that one may calculate the 
Chern number of the occupied quasi-energy bands, as for 
the unperturbed system. This works if there are clear 
gaps between the quasi-energy bands. In the context 
of this work, this will work in general at high frquen- 
cies. In most other cases, even though the analysis in 
the stripe geometry reveals edge states, these are often 
mixed with continuum states that close the gaps between 
bands, preventing the calculation of the Chern number 
of these bands. Also, due to the periodicity in Floquet 
quasi-energy space, we may have bands with Chern num¬ 
ber that vanishes if the number of edge modes that enter 
the band equals the number of modes that leave ifP. 


( 'u(/ta,, 1, u(^kxj I5 '!') u( kxt I5 d) u( kx^ f 5 '!') ‘ ‘ ^ 


The time evolved states due to the perturbation are the 
result of the diagonalization of the truncated Floquet ma¬ 
trix and each component {U = u or U = v) is of the form 

l^{kx,jy,(J,t) = '^e^'^‘^*(j)raikxJv,cr) (17) 

m 

The charge current operator along direction x at a 
given position jy along y is given by 

) (n ^ Vi/d . ( -^'sin(fc,) -^acos{kx) \ , 
3c\3yj ^ 2^ k^,jy ^acos{kx) —isin{kx) ) 

(18) 

where'0^ ^ ^ j • The current has con¬ 

tributions from the hopping and the spin-orbit terms. 

One may also define a longitudinal spin current, js(,3y), 
taking the difference between the two diagonal compo¬ 
nents of the charge current. The other terms correspond 
to spin-flip terms and do not contribute to the z compo¬ 
nent of the spin current. 

The average value of the charge current in the ground- 
state is given by summing over the single particle occu¬ 


Since our aim here are the currents associated with the 
edge modes, it is more enlightening to look at the band 
structure in the stripe geometry, since both the edge and 
the bulk states (projected along one spatial direction) are 
found. 


V. CURRENTS 

We consider a finite system of dimensions Nx x Ny. We 
apply periodic boundary conditions along the x direction 
and use a momentum representation and open boundary 
conditions along the transverse direction, y, solving the 
problem in a ribbon of width Ny. Writing 

'4’k^,jy,a = ~/^ ® ’ 

Jx 

we may rewrite the Hamiltonian matrix in terms of these 
operators. The diagonalization of this Hamiltonian in¬ 
volves the solution of a {4:Ny) x (iNy) eigenvalue problem 
for each momentum kx ■ The energy states include states 
in the bulk and states along the edges. The eigenstates 
give us directly the wave functions in real space. 

In this representation the states are column vectors 
that with no external time dependent perturbation is of 
the type 


’x^A^yiT) '^{kx •, Ny ^ i') v( kxT Ny^ 'I') v( kx, Nyj .1.) ^ 

,_ m 


pied states (negative energies) in the usual way 

3c{jy) = {jc{jy)) — ^ ^ 
kx ,71 

sin kx ['l^n ( jy-} T)^n ( - jy ; T) 

“1“ '^n( ; Jy;( ^x-jyi 

ItOi 

^ COS kx 1^71 ( ^x j jy 5 T )^n { -jy : 4') 

'Cn( kx^jyjl-)v^{ ^T-JyiT)]} (19) 

Here the functions are of the type 

Un{kx,jy,Cr) = jy,Cr) ( 20 ) 

m 

where as usual a =t, i- 

It is also interesting to consider the magnetization of 
each momentum value. The spin polarization is obtain¬ 
able by 

m{kx,3y) = (21) 
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and the average spin polarization in the groundstate is 
given by 

m{kx,jy) = 

n 

{~in{~kx, jin t)^n{~kx-jy, t ) 

+ Vn{ kx, jy, kx-jy,i-)} 

( 22 ) 

A. Charge currents 

The charge currents at the edges of the unperturbed 
Hamiltonian are well understood. If there is TRS the 
currents vanish and if TRS is broken the edge charge 
currents are finite (finite Chern number). We consider 
as examples a dx,dy triplet superconductor and a, p + ip 
triplet superconductor. In the first case, if the magnetiza¬ 
tion vanishes, the system has TRS and vanishing charge 
edge currents in the topologically trivial phases. The 
charge current also vanishes in the Z 2 topological phase 
but the spin edge currents are non-vanishing. In the case 
of the p + ip triplet superconductor there is no TRS and 
the charge edge currents are finite. 

We consider in the case of the dx,dy pairing a set of 
parameters d = 0 . 6 , = 0 . 1 , = 0 , a = 0.6 and dif¬ 

ferent values for the chemical potential and the magne¬ 
tization. In the case of the p + ip pairing we consider 
d = 0 , As = 0 . 1 , dz = 0 . 6 , a = 0 . 1 . 

In Fig. I^we show results for the profile of the charge 
current as a function of y for the various cases. We com¬ 
pare the unperturbed case with the perturbed one by 
considering that at w = 6 it is enough to truncate the 
Hamiltonian matrix at / = 2. To calculate the currents 
we sum over the states in the first Floquet zone. Also 
the results are for time t = 0 or any multiple of the time 
period T. 

The results for the dx, dy pairing consider four points in 
the phase diagram: namely a trivial phase /r = — 5, = 

0 (case Pi with C = Q), a Z 2 phase p = —3, = 0 (case 

P 2 with (7 = 0), and two points in Z topological phases 
P 3 with p = —5,Mz = 2 (with (7=1) and P 4 with 
p = —l,Mz = 2 (with (7 = —2). The results for the 
p + ip pairing consider P 5 with p = —5, Mz = 0 (with 
C = 0), Pq with /i — 2, Mz = 0 (with C = 2), Pj with p = 
—3,Mz = 2 (with (7 = 1 ) and Pg with p = —l,Mz = 2 
(with (7 = —2). These points are indicated in Fig. 

In Fig. we consider a periodic driving in the chem¬ 
ical potential with amplitude pd = 1 and in Fig. 
in the magnetization with the same amplitude. Focus¬ 
ing first on the unperturbed case we see that the current 
at the border vanishes for Pi (C = 0) and that for P 4 
with (7 = —2 the current is larger at the edge, but de¬ 
cays faster with distance from the border with respect 
to the P 3 , (7 = 1 case. Adding the perturbations in the 
case of the trivial Pi case, there is no charge current if 
driving the chemical potential but if one drives the mag¬ 
netization there is a small current close to the edge. The 




' ^ ^ ^ 

1 
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FIG. 6 : (Color online) Charge current profile for I — 2 for 
a dx,dy superconductor for Pi,P 3 ,P 4 for (a) pd = 1 and (b) 
Mzd = 1. Only one half of the system is shown since the cur¬ 
rent profile is anti-symmetric around the middle point. The 
system size is Ny = 100. Charge current prohle for I = 2 for a 
dz superconductor for P 5 , Pg, P 7 , Pg for (c) pd = 1- Note that 
the current profile has no symmetry around the middle point 
in some cases. The results for the dx,dy pairing consider three 
points in the phase diagram: p = —5, Mz = 0 (case Pi with 
C = 0), P 3 with p = —5, Mz = 2 (with (7=1) and P 4 with 
p = —1, Mz = 2 (with C = —2). The results for the p + ip 
pairing consider P 5 with p = —5,Mz = 0 (with C = 0), Pg 
with p — 2, Mz = 0 (with (7 = 2), P 7 with p = —3, Mz = 2 
(with (7 = 1) and Pg with p = —l,Mz = 2 (with (7 = — 2 ). 
These points are indicated in Fig. 
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FIG. 7: (Color online) Influence of the spin-orbit coupling on 
the spin current in the unperturbed case. The parameters are 
d — 0.6, As = 0.1, dz = 0, Mz = 0, /r = —3, —5. Increasing the 
spin-orbit in the topologically non-trivial phase is detrimental 
while in the topologically trivial phase, increasing the spin- 
orbit coupling considerably increases the spin current. 

effect of the perturbations on the topological cases is how¬ 
ever quite significative, particularly in the case of C = 1, 
with a current inversion with respect to the unperturbed 
case. A similar result has been found in the context of 
dx 2 -y 2 + idxy pairinj^. In both cases the spatial extent 
of the currents away from the border is considerably in¬ 
creased. However, we have checked that integrating the 
currents over y from the edge to the middle point the to¬ 
tal current is approximately the same. Also, we see that 
if the unperturbed state is topologically non-trivial the 
effect of the two types of perturbations is very similar. 

In Fig. we show results for the p + ip pairing. Ex¬ 
cept for the case with (7 = 0, the charge currents are 
finite in the unperturbed and perturbed cases. We note 
that, even for the (7 = 0 case, there is a small charge 
current at the border, probably due to some finite size 
effects. In general, the driven chemical potential consid¬ 
ered here increases the charge current. In the cases with 
zero magnetization the current decreases from the bor¬ 
ders but at finite Mz the current is finite throughout the 
whole system and is not symmetric around the middle 
point. 

B. Spin currents 

In this subsection we calculate the longitudinal spin 
currents at the edges of the system, both for the unper¬ 
turbed and perturbed systems. 

Consider first the unperturbed triplet superconductor. 
We start with the dx , dy triplet superconductor in a triv¬ 
ial phase in zero magnetic field. In this regime there are 
no edge states and the system has a finite gap around zero 
energy. In Fig. [^we consider the influence of the spin- 
orbit coupling on the spin current in this trivial phase 


FIG. 8: Influence of intensity of periodic driving for the cases 
of p,d,ctd, Mzd on the spin current (Z = 2). The parameters 
are p, = —5,d — 0.6, As = 0.1, dz = 0,Mz = 0,a = 0.6 and 
the frequency of the periodic driving is w = 6. Note the re¬ 
versal of the direction of the spin current at small couplings. 
At weak coupling the increase of the spin current is linear. As 
the transition from weak to strong coupling occurs the spin 
current becomes non-linear. Changing the spin-orbit coupling 
has a small effect on the edge spin current while both chang¬ 
ing the chemical potential and the magnetization have strong 
effects includind a change of the current direction. 


and compare it with the case of a topological Z 2 phase, 
where there are edge states that, due to the spin lock¬ 
ing, counterpropagate at each edge yielding a null charge 
current but a finite spin current. This is shown in the 
results for the case with p = —3, Mz = 0. Increasing the 
spin-orbit coupling decreases the longitudinal spin cur¬ 
rent. An opposite effect is observed in the trivial phase 
for which, as the spin-orbit coupling increases, the spin 
current increases. At very small coupling the spin current 
is very small but non-vanishing, probably due to finite 
size effects. Even though there are no edge states the 
spin current is also carried out by the states in the con¬ 
tinuum above and below the gap. When there are edge 
states these carry most of the current, but there is also 
a contribution from the other states in the continuum. 

Consider now the perturbed cases. There are two 
parameters that characterize the perturbation: the fre¬ 
quency of the periodic driving and its amplitude. As 
discussed above, small frequencies require a large matrix 
to be diagonalized, particularly to determine the edge 
states. Very large frequencies are amenable to analyt¬ 
ical solutions using the Magnus expansion and can be 
interpreted as a renormalization of the parameters of the 
original Hamiltonian, or by the addition of new terms 
allowed in the original unperturbed Hamiltonian. The 
cases treated here respect to intermediate values of the 
frequency, to allow for a simple truncation scheme and, 
at the same time, a resonance condition between bands 
that is most interesting in changing the topology of the 
system. Most results reported here are for a small to 
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FIG. 9: (Color online) Spin current for the triplet superconductors dx,dy ((a), (b)) and d^ ((c), (d)) for the periodic drivings 
/id = 1 ((a), (c)) and = 1 ((b), (d)). 
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FIG. 10: Time evolution of the spin current for a dx,dy triplet superconductor for the periodic drivings (a) ^d, (b) ad, and 
(c) = 1 for Pi, P 2 , Pi, Pa, respectively. 


moderate amplitude. If the perturbation has a very small 
amplitude, it can be dealt with using perturbation the¬ 
ory. Very large amplitudes deviate considerably the sys¬ 
tem from its initial state. Therefore we have considered 
here moderate amplitudes. To illustrate the effect of the 
increase of the perturbation amplitude we show in Fig. 
for the same triplet superconductor, the effect of chang¬ 


ing the amplitude of a perturbation in the chemical po¬ 
tential, spin-orbit coupling and magnetization. As above 
for the charge currents, the results hold for time t = 0 or 
for any multiple of the period T. The results for the spin 
current are obtained integrating the local spin current 
accross half the sample. 

In all cases there is a sign reversal of the spin current. 
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At small couplings, the behavior of the spin current {I = 
2 ) is approximately linear and the spin current does not 
change appreciably. Increasing the amplitude in the case 
of the spin orbit coupling, ad, does not affect much the 
spin current even though the response is no longer linear. 
On the other hand, the change in the chemical potential 
leads to a significant change of the spin current, changing 
back its direction for large enough values of the chemical 
potential, fj,d, and the effect is even more pronounced in 
the case of the magnetization, M^d- These results were 
obtained taking the unperturbed as a trivial phase with 
a very small spin current. The spin current is therefore 
significantly enhanced by the periodic driving. 

In Fig. [^we show results for t = 0 (or a multiple of the 
period) of the spatial profile of the spin current for the 
triplet superconductors dx,dy ((a) and (b)) andp+ip ((c) 
and (d)) obtained changing the chemical potential, pd, or 
the magnetization, Mzd- In the case of the first pairing 
symmetry, the spin current is antisymmetric around the 
middle point and so it is enough to consider half of the 
system. In the case of the p + ip superconductor the 
current has no symmetry and so the full profile is shown, 
as before for the charge current. 

The results are to some extent qualitatively the same 
as for the charge current, in the sense that the periodic 
driving induces a large spin current, if the unperturbed 
system is trivial, and enhances the spin current if the sys¬ 
tem is originally topological. The spin current oscillates, 
there is also reversal of direction due to the periodic driv¬ 
ing, in the case of C = —2, and the spatial extent of the 
current profile increases. As for the unperturbed case, 
when there are edge states, the largest contribution to 
the current originates in them, but the continuum states 
also contribute. 

In the case of the p + ip superconductor the currents 
are spread along the whole system and in general are not 
symmetric around the middle point of the y direction. 

At arbitrary time values the currents oscillate. In Fig. 
| 10 | we show the time evolution of the spin current within 
one period, T, for the case of the dx,dy triplet super¬ 
conductor driven by the chemical potential, spin-orbit 
and magnetization. We consider as unperturbed states 
in each panel the trivial case Pi , the Z 2 topological point 
P 2 , and the Z topological phases P 3 and P 4 . In the triv¬ 
ial phase, at short times, the perturbations decrease the 
spin current with respect to the initial time value (or any 
multiple of the time period). In the non-trivial phases 
the chemical potential and magnetization perturbations 
increase the spin currents at short times with respect to 
the initial time value while the spin-orbit coupling acts 
the opposite way. 


C. Spin polarization 

Associated with the spin currents in the topological 
phases, it has been shown that the eigenstates have non 
trivial spin polarizations that depend strongly on the 



FIG. 11: (Color online) Spin polarization of unperturbed 
(Z = 0) and perturbed (Z = 2) system with /id = 1 for (a) and 
(b) Pi (trivial case) and (c) and (d) P 2 (Z 2 case). In these 
cases C = 0 and Mz = 0. The maximal values of the spin 
polarizations are (a) 7.3 x 10~®, (b) 1.9 x 10~®, (c) 2.5 x 10“^ 
and (d) 5.8 x 10“^. The maximal value displayed in panel (c) 
has been reduced to 0.007 for better visualization. 


momentunPSI. It was shown that, particularly for cases 
where one has flatbands, there is a strong polarization 
effect. Here we will look at the spin polarization of the 
Floquet states and compare with the unperturbed cases. 

In Figs. and we show the spin polarization as 
a function of momentum for each quasi-energy state for 
the various examples, Pi, P2, P3, P4 of the dx,dy super¬ 
conductor, both for the unperturbed and the perturbed 
cases with I = 2 . 

The unperturbed cases in Fig. El show the nontriv¬ 
ial momentum distribution of the spin polarization for 
cases where (7 = 0 and Mz = 0. Summing over all mo¬ 
menta and states, the total polarization vanishes. The 
spin polarization shown in Fig. [TT}r has higher values at 
the edges of the gap and opposite signs as the momen¬ 
tum changes sign, due to the time reversal symmetry. 
Note, however, that the spin polarization is very small 
(^ 10“^). In this regime of parameters the system is in 
a trivial phase. 

Adding the periodic driving edge states are generated, 
as shown above, the spin polarization has higher values 
that are more uniformly spread over the various energy 
states and momenta, and with clearly higher values at 
the induced edge states at the border of the Floquet zone, 
again with opposite signs as the momentum changes sign 
since time reversal symmetry is not broken. However, the 
spin polarizations have a complex structure over the con¬ 
tinuum. Panels (c) and (d) refer to a phase that shows 
edge states already at the level of the unperturbed Hamil¬ 
tonian but with time reversal symmetry. The spin polar¬ 
ization is much stronger and sharply localized along the 
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FIG. 12: (Color online) Spin polarization of unperturbed 
(Z = 0) and perturbed (Z = 2) system with = 1 for (a) 
and (b) P 3 {C = 1) and (c) and (d) P 4 (C = —2). In these 
cases G 7 ^ 0 and Mz 7 ^ 0. The maximal values of the spin 
polarizations are (a) 3.8 x 10~^, (b) 7.8 x 10~^, (c) 4.7 x 10“^ 
and (d) 8.3 x 10“®. The maximal values of the panels have 
been reduced to 0.005. 


states that constitute the edge modes. Turning on the 
perturbation, the weight along the original edge states 
remains evident, but is reduced and there is a pile up of 
spin polarization along the new Majoranas generated at 
the edge of the Floquet zone. Note that, even though 
there is an alternancy of the polarization sign when the 
momentum sign is reversed, the two edge states are quite 
visible. 

The results in Fig. 12 are different. They apply to 
cases for which the Chern number does not vanish, TRS 
is broken in the unperturbed system and there is a finite 
magnetization. The spin polarization is once again con¬ 
centrated mostly along the edges but only one edge state 
is visible. The edge may be switched to the other one by 
reversing the sign of the magnetization. If the edge states 
are present in the unperturbed case, it is also clear that 
adding the time perturbaton decreases the polarization 
overall amplitude. Considering unperturbed topological 
phases, the spin polarization is highly concentrated on 
the edge states and less so in the perturbed case. In 
any case there is an enhancement with respect to the 
bulk states along the edge states, with either positive or 
negative spin polarizations in both branches of the edge 
states. However, in the presence of a magnetic field this 
symmetry is gone. 

In the case of the p + ip superconductor (not shown 
here) the states of the unperturbed Hamiltonian below 
the gap have a spin polarization that is large (and posi¬ 
tive) for the negative momenta values but is considerably 
smaller for the positive momenta. The opposite happens 
for the positive energy states. Both cases correspond to 


zero magnetization but the chiral nature is evident, both 
in the trivial case (7 = 0 (expressed in the bulk states) 
and in the topologically non-trivial phase {C = —2), (also 
in the edge states). This effect is somewhat washed out 
turning on the periodic driving, but both for the unper¬ 
turbed and perturbed cases, only one edge state has a 
large spin polarization. In the perturbed case the edge 
state is immersed in the continuum. Nevertheless the 
spin polarization clearly singles it out. 


VI. CONCLUSIONS 

A periodic driving on a topologically trivial system in¬ 
duces edge modes and topological properties. In this 
work we considered triplet and singlet superconductors 
subject to periodic variations of the chemical potential, 
spin-orbit coupling and magnetization in both topologi¬ 
cally trivial and nontrivial phases and studied their in¬ 
fluence on the charge and spin currents that propagate 
along the edges of the two-dimensional system. The most 
relevant case is the generation of the charge and spin edge 
currents due to the periodic driving in an otherwise triv¬ 
ial system. In some cases the generated currents are quite 
high and higher than in the unperturbed case providing 
a way to generate and control high spin currents with po¬ 
tential interest in the context of spintronics, particularly 
since some of the driving protocols involve electric or gate 
potential means, instead of a direct magnetic manipula¬ 
tion. Starting from a topologically non-trivial phase, the 
edge states originally present at zero energy are comple¬ 
mented by edge states at the edge of the Floquet zone at 
finite quasi-energies close to ±ry/2. In general, the peri¬ 
odic driving smears the edge states, as shown particularly 
in the study of the spin polarization of the quasi-energy 
states, with a complex structure both in the unperturbed 
and perturbed cases. The chirality of the edge states is 
particularly seen when TRS is broken either by the pres¬ 
ence of a finite magnetization or due to the intrinsic TRS 
breaking in a. p ip triplet superconductor. 

The transport signatures of edge states in driven 
topological insulators or due to the Floquet Majorana 
fermions in driven topological superconductors, involve 
the presence of leads to inject and collect charge or spin 
currents and therefore require the coupling of the sys¬ 
tem to outside reservoirs, with the associated issue of 
dissipatioiP^. A calculation of the differential conduc¬ 
tance has found a quantized conductance sum rul^^, 
which generalizes the quantized zero-bias conductance 
contribution of the Majorana fermions when the system is 
not periodically driven. Here we neglected any couplings 
to external reservoirs. The charge currents calculated in 
this work may be detected considering the magnetic field 
they generate in their vicinity, as proposed in the context 
of unperturbed systems in Ref.^. Barring in mind the 
difficulty in detecting experimentally the edge currents, 
it may be challenging as well to detect the generated cur¬ 
rents considered here. 
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